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. Abstract. We study the timc evolution of the quantum-classical correspondence 

(QCC) for the weU known model of quantised pertm'bed cat maps on the torus in 
the very specific regime of semi-classically small perturbations. The quality of the 
QCC is measured by the overlap of classical phase-space density and corresponding 
Wigner function of the quantum system called quantum-classical fidelity (QCF). In 
55 ■ the analysed regime the QCF strongly deviates from the known general behaviour 

^ [ discussed in [1], in particular it decays faster then exponential. Here we study and 

qh' explain the observed behavior of the QCF and the apparent violation of the QCC 

principle. 
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The quantum-classical correspondence (QCC) is tlie basic principle underlying any 
physical quantisation of a classical system. According to this principle the quantum 
system should behave similar to the corresponding classical system with increasing 
5^ i energy or decreasing effective Planck constant. The importance of QCC as a tool in the 
study of quantum systems was recognised very early in the development of quantum 
mechanics with the Ehrenfest theorem and later by the introduction of semi-classical 
methods [2j. The study of QCC gave in the 80's birth to quantum chaology - research 
area devoted to study the connections between dynamical properties of classical systems 
and corresponding quantum systems 

The QCC can be explored and discussed using various tools and methods available 
in the theory of classical/quantum systems. In a recent paper [T], a phase space 
representation has been used to study the time evolution of QCC in generic chaotic 
systems on compact classical phase space. 

The QCC is there quantified using the so called quantum-classical fidehty (QCF), 
namely the integrated overlap between the classical phase-space density and the 
corresponding Wigner function. It has been shown that in classically chaotic systems, 
after some initial plateau, the QCF decays exponentially in time with decay rate 
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coinciding with the maximal Lyapunov exponent A. While it is common knowledge 
in the realm of quantum chaos that the phase space correspondence between classical 
and quantum mechanics drops down on the scale of Ehrenfest time t_B ~ —\ogh/X, 
the exact dependence of initial plateau on dynamical properties and on Hilbert space 
dimension are still important open questions that we aim to address here, at least 
in a special case. In particular, here we discuss QCC using QCF for the so called 
perturbed Arnold cat map[lj on a torus = [0, 1]^. The (unperturbed) cat map is a 
paradigmatic example of a classical uniformly hyperbolic chaotic systems. It has been 
one of the first extensively studied quantum maps and since then it has been used 
several times to prove or disprove various conjectures concerning statistical properties 
of eigenfunctions for quantum system with strongly chaotic classical motion (see for 
example |6] and references therein). Because of the linearity of the classical motion, 
its quantum counterpart inherits a natural non generic number-theoretical structure, 
reflected for example in the rigid distribution of eigenvalues and also in the so called 
ezactness of the Egorov property, which roughly means that classical and quantum time 
evolution perfectly commute. A generic behaviour of eigenvalues, i. e. a good agreement 
with the predictions of random matrix theory, can be gained by perturbing the linear 
cat dynamics by composing the map for example with a time one flow generated by a 
given (global) Hamiltonian over the torus [71 [8] . 

We are here interested in exploring the time evolution properties of these perturbed 
maps. In particular we aim to study how QCC decays in the presence of perturbation 
and especially the time scales of the initial plateau in QCF. 



2. Quantisation on a torus and Egorov property 

We recall here the basic facts of quantum mechanics over the torus which we need 
in the paper, see [6] and [9] for further details. The system is quantised on a torus 
by introducing a position basis {\qn) '■ Qn = -^}neZN a momentum basis 

{\Pm) '■ Pm = ^}meZ]v Hilbcrt space Hn dimension A^. We apply periodic 

boundary conditions \qn+N) = kn) and \pm+N) = \Pm)- The two basis are related by 
the discrete Fourier transform: 

(g„b„) = -^e^t-. (1) 

Then according to Weyl-Wigner quantisation we associate an operator A to a classical 
observable a defined over grid points = {xn,m = (2^; i^)}(n,m)Gz2j^ classical 
phase-space using following relations 

(^n,m ~ tr{y4.yl.„ „jj- , A = Qv/{0') — ^ ^ 0,n,m-^n,m j (2) 

where An,m is called the point operator or the kernel of the Weyl-Wigner formalism 

= — = Yl e~''^'n<ln-k){qk\ , (3) 
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We say that a is phase space representation of operator A or A is quantisation of phase- 
space function a. The phase-space representation of a density operator corresponding 
to a pure state p = \ip){4'\ is the Wigner function W^{n,m) defined as 

W4n,m) = —j= V (^|g„-fc)(g;^|^)e-t'=- , (n, m) e Z^^ . (4) 

with normalisation J2{nm)€^^ W^^n,!!!,)^ = 1. Namely, 

Qw(W^^)=P=|^)(^|. (5) 

Let us assume that M : — is a classical discrete, area preserving, map on the 
torus. Then it is possible to associate to M a corresponding quantum evolution operator 
U : TLn ^ 'Hn- In the foUowing, is always an even integer in order to avoid certain 
technicahties in quantisation (see [6] and references therein). The propagator U will 
satisfy an Egorov estimate, namely [TO] : 

hm II f/tg„(a)f/-Q„(aoM) 11=0. (6) 

A'"— >oo 

In particular, if the classical map is a linear automorphisms i. e M(x) = M.x, where 
matrix M G SL{2, Z) (e. g. the cat map) then Egorov is exact 

U^Q^{a)U = Q^iaoM). (7) 
3. Cat map and breaking of Egorov property 

The classical dynamics over the torus that we study here is given by the map M : — »• 

I V ^ p' = P + kq + e-Y{q) mod 1 

g,p = M g,p , , (8) 

q = q + p moa 1 

with A; G N, perturbation parameters e = (eo, ei, £2) G M.^ and perturbation function 

V(g)= Qg',-^cos(27rg),g^ . (9) 

The quantum evolution operator corresponding to the perturbed cat map U can be 
written as 

U = exp (-i^^^) exp (i^^^^ + iA^e ■ , V = V (^^^^ . (10) 
where we have for convenience introduced auxiliary operators 

n\qn) = n\qn) , m\Pm) = m\pm) ■ (11) 

If the perturbations are neglected e = O, we obtain the usual linear cat map system, 
where the map and evolution operator are denoted by 

M, = M\e=o, U, = U\e=o- (12) 
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As already remarked, it is well known that this system is Egorov exact. Moreover, the 
classical cat-map is uniformly hyperbolic with the Lyapunov exponent 



A(A;) = log 



(13) 



For the coming analysis, it is convenient to write the classical map and the quantum 
evolution operators of the perturbed cat map as 

U = Uc exp (^iNe ■ v) , M = Mc + (e ■ V, e • V) . (14) 

We compare the classical and quantum evolution of the perturbed cat map in the 
classical phase space at some fixed dimension N and perturbation e. The e = ||e|| 
is refered to as perturbation strength. Here we are mainly interested in the particular 
case of semi classical small perturbations A^e ^ 1. More precisely, the classical system 
starts from a smooth probability distribution p : ^ M resembling a Gaussian packet 
on phase space at the point {qo,Po), 

P(,,,„)(g,p) = D^,(,,„„) [j^e--^^--^^)^) [j^e--^^--^^)^) . (15) 

where the scalar factor is pinned down by the normalisation 

Yl P(go,Po)K™) = l' (16) 



(n,m)GZ? 



'2JV 



which has a simple leading term in the asymptotic approximation, N oo, reading 

X [(1 + 2 cos(47rArgo)e--^) (l + 2 cos(47riVpo)e-"^)] . (17) 



The quantum counterpart is initially in a coherent state |0) with a Wigner function 
similar to the classical distribution (see [Tl [TT]): 

W^n, m) = p{xn,m) + e-\^^^^\ , (n, m) G Z^^ . (18) 

We then let these two systems evolve up to time t G Z* using equations 

p* = p o , 10*) = f/*|0) , p* = U'pU-' , (19) 

and observe QCC between these two systems by calculating the overlap of the density p* 
and corresponding Wigner function W^t . The overlap is called quantum-classical fidelity 
(QCF) defined as 

m= E W^.{n,m)p\xr.,J=tT{p'Q{p')] < 1 + e-'^^^^l (20) 

Because the perturbed system is not Egorov exact, the QCF decreases with time. In 
figure[T]we show numerically obtained the decay of average QCF {F{t)) for different k 
and perturbation strengths e using perturbation vector e = (e, 0,0), where (•) denotes 
the uniform average taken over the initial positions of the coherent packet. The QCF 
does not decay up to some time called the breaking time tbr, which increases with 
decreasing perturbation. Beyond tbr the QCF decays "very fast" (as we will argue. 
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Figure 1. The average QCF F{t) with perturbation vector e = (e, O, 0) for different 
perturbations strength e at Hilbert space dimension N = 512 and k — 1,2 (a, b). The 
average is taken over 100 initial Gaussian packets uniformly scattered over phase-space. 



faster then exponential) and eventually converges to the ergodic plateau given by 
the decay is in fact visually faster than the generally expected exponential Lyapunov 
decay {F{t)) ~ exp(— At) [1], which is inserted in the figures. Basically the same scenario 
happens in other choices of perturbation vectors e as we can see in figure [2j In the case 
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Figure 2. The average QCF F{t) for different perturbation strengths e with 
perturbation vector e = (0,e,0), (0,0, e) (a, b) at Hilbert space dimension N = 512 
and k — 1. For averaging see caption of fig. [TJ 



of constant classical perturbation e = (O, O, e) shown in figure [2lb, the correspondence is 
broken mainly by a rigid shifting of the deformed packets in the quantum and classical 
picture. Therefore the convergence towards the ergodic plateau is less smooth as in 
other cases. At this this point it is difficult to deduce the correct functional form of 
the QCF F{t). Nevertheless, in the foUowing we present a theoretical explanation of 
these numerical observations, disclosing the super-exponential nature of the QCF decay 
in this particular regime of perturbation. 
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We are interested in the evolution of the QCF in the limit of small perturbations 
e — > 0. In this regime we examine the time tbr(p) on which average QCF {F{t)) drops 
below some value p: 

tM=mm{teZ* -.{Fit)) <p}, (21) 

where the average (•) is taken uniformly over positions of the initial Gaussian packets. 
It is meaningful to express the dynamics relative to the unperturbed cat map writing 

p* = p* + 6p' , p' = pi + Sp' , (22) 
where the dynamics of the cat map case is given by 

pI = U'JU-' , p* = p o . (23) 
and due to Egorov property these are connected by 

Pc = Q(Pc) • (24) 
By inserting ansatze fl22l) into the formula fl20l) we obtain QCF expressed in terms of 
deviations from the unperturbed case: 

F{t) = 1 + tr {6p'Q{pi)] + tr {plQ{6p')} + tr {6p'Q{6p')} . (25) 

Due to existence of the Egorov property in the cat map, the approximated QCF can be 
expressed in terms of the quantum fidelity Fq(t) [12] and the classical fidelity Fdt) 



as 



F(t) = |Fq(t)|2 + Fc(t)-l + tr|5p*g(5p*)} , (26) 



where F^ and F^ are here written as 



F^{t) = m-'tj'M , ^c(t) = E p(M-*(x.,^))p(M-*(x.,^)) .(27) 



The relation 0261) is very instructive and helps to understand the behaviour around 
the initial plateau, but it seems to us that the study of the plateau itself was greatly 
avoided in the past. In the foUowing we discuss the second term and the third term in 
( 125|) denoted by 

/l = tr|p*Q(5p*)} = Pl{^n,m)5p\xn,m) , (28) 

(n,m)GZ2jY 

/2= tr{5p*g(p*)}=tr{5pV*} . (29) 

The last term in fl2Sl) and fl2B]) are the second order corrections, which we do not 
discuss in detail. In order to understand /i ( l28l) we discuss the deviation between 
trajectories of a chaotic and ergodic map = M^^ : — and of its perturbation 
+ 50 = (M + 5M)~^, starting at the same point x. The deviation is defined as 

50t(x) := (0 + 6<P)\x) - (t>\x) , 0*+i(x) = 0*(0(x)) . (30) 

and obeys in the limit 50 O the following recursion 

5(t>t+i{x) = (0 + 50)(0*(a;) + 5ct>t{x)) - 0*+i(a;) , (31) 

= (V0)(0*(x))50,(x) +50(0*(x)) , (32) 
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where we have neglected second order corrections. By iterating this equations from a 
given initial position x, the deviation is written as a series 



6u^) = Y, 



k=0 



n(v0)(0'(x)) 



.l=k 



50(0^-1(0;)) + (33) 



Then by taking into account that map is chaotic and ergodic with Lyapunov exponent A, 
we get in the limits t — oo and 6(f) O, applied in given order, the leading contribution 
of the deviation expressed as 

(50i(x) = 0(50) e"* , (P0t(x)l|) ^ eAe^' . (34) 

where (•) denotes the uniform average over initial positions x. The constant A G M 
depends only on the type of perturbation and dynamical properties of the map. By 
plugging this result into expression /i (1281) we obtain 

h= ^ P(a;„,m)P {xn,m + S4>ti^n,m)^ " 1 , S4>t = ^<Pt ° M* , (35) 

where we have used that the cat map M conserves the grid Gn'- M{Gn) = Gn- By 
taking into account the exphcit form of p f|T5|) and considering only the behaviour about 
the central point of the Gaussian packet the above expression is approximated as 

h ^ exp (-7riV||50,(go,Po)f ) - 1 , {h) ~ exp(-7riVA2e2 exp(2At)) , (36) 

with (•) representing the uniform average over position of the initial coherent packet. 
The approximation is meaningful up to times eA^^ exp(At) = 0(1), when deformation 
of the packets can be neglected. This is especially appropriate to described the case of 
constant classical perturbation. In the limit of small perturbations the leading term in 
expression /i scales with time and perturbation as 0(A^e^ exp(2At)), where the changes 
of QCF are small. The behaviour of the expression /i is obtained by considering the 
fact that 

t 

U' = lJl + iNe-Ylj'^\Ul-^ + 0(^{N\\eV\\ft) , (37) 



fc=i 



which yields 



5 p' = iNe . [s'pI - p-'S-'] + O [{N\\eV\\)H) , (38) 
t 

S' =Y1 ^c"' Vf/,^ . (39) 



k=l 

By plugging this into I2 ( 129|) we get 



h = 2Ne.Q{tT{s'pi}}+0[{N\\eV\\?t) • (40) 

By assuming that lim^^oo t~^S' 7^ O we see that the leading term in /2 scales as 0{Net) 
in time. Then by considering results /1 (!36|l and /2 ( 140|) we get the leading order 
contributions to QCF reading 

F{t) ^2Ne-'^ |tr {s*p*} } + exp [-7rN\\6UQo,PoW) • (41) 
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In the limit of the small perturbations the last term in fl41l) is the dominant. This 
is supported also numerically as we demonstrate in the figure [3], where we show 
G(t) = log(— log((F(t))) for different perturbation vectors e, perturbation strength e 
and k. We see that the average QCF evolves foUowing the curve G(t) ^ 2\{k)t + const. 
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Figure 3. The representation of average QCF {F{t)) evolution calculated using 
perturbation vector e equal to (e, 0,0) (a), (O, e, 0) (b) and (0,0, e) (c) for different 
e and /c = 1, 2 at the dimension N = 512, where data presented in figs. [T]and[2]is also 
considered. For averaging see caption of fig. [TJ 

according to the dominant term in QCF Ii{t) (136!) almost up to the time, when QCF 
intersects the ergodic plateau given by G{t) k, log log A^. The plot G{t) has an initial 
plateau due to finite arithmetic. We conclude that the QCF decays in average towards 
the ergodic plateau faster than exponentially as 

(F(t))=exp(-|0(exp(|0(t)|))|). 

The expression for QCF fj4T|) obtains in the limit 
form 

(F(t)) = 1 - 0{Ne^e^^') + 0{Net) . 

where the first non-constant term is dominant in F {t). In this perturbation approach 
we can approximate tbr (p) (121]) for fixed 1 — p ^ 1 as 
log(l - p) - log(iVe2) 



(42) 

A^eexp(At) ^ 1 a simple scaling 

(43) 



t 



2A 



(44) 
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which in the limit of infinitesimal perturbations obtains foUowing asymptotic form 

AtbrX-loge, e^O. (45) 

We see that at fixed and p the time depends only on Lyapunov exponent A and 
perturbation strength e. 

4. Numerical result on the breaking time 

In the foUowing we present numerical results of the breaking time tbr in our perturbed 
cat map. We explore in particular its dependence on the perturbation strength e and 
the Hilbert space dimension A^. 

The figures H] and [5] show plots of tbr in dependence of e for all three types of the 
perturbations. Because we are discussing a discrete dynamical system, the break time 
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Figure 4. The dependence of ^br on perturbation strength e in the cases fc = 1 (a) 
and fc = 2 (b) by using perturbation vector e = (e, O, 0) at = 512 . 



thr{^,p) is a discrete function of e G M^. In figure H] we show tbr as a function of e in the 
case of non-smooth perturbation e = (e, O, 0) for two values of the classical parameter 
k. In order to improve representation we show plots for several p at the same time. We 
see that the heuristically obtain formula Atbr ~ ~" log ^ fit perfectly onto the numerical 
results. The dependence of tbr on e in the presence of smooth perturbations e = (O, e, 0) 
and e = (O, O, e) is shown in figure [51 We notice that the gross dependence of the break 
time is basically independent of perturbation. 

The break time tbr ( I2T1) depends also on the Hilbert space dimension N. In the 
limit of small perturbations eN <C 1 we obtain from (jH]) the foUowing dependence on 
N: 

Atbr ~ const. — ^ log , (46) 

where the constant depends on p, A and details of the initial packets. The numerical 
results shown in figures [6] and [7] in the case of using smooth and non-smooth 
perturbation, respectively, confirm the theoretical dependence. But due insufficient 
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Figure 5. The dependence of tbr on perturbation strength e in the case of pertbation 
vectors e = (O, e, 0), (O, O, e) (a,b) at = 512 and k = 1. 
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Figure 6. The dependence of ibr of Hilbert space dimension N at perturbation 

e = (10^1°, 0, 0) in the case fc = 1 (a) and fc = 2 (b). 



range in variable logA^ we can not check the prefactor in scaling relation (H6l) very 
accurately. 

Notice that the break time tbr is decreasing with increasing A^. At the first look 
this would appear as a contradiction to the known QCC principle, which states that the 
quantum system should behave as classical system in the hmit N ^ oo. But this is not 
the case: with increasing N eventually Ne ~ 1 and the perturbation approach becomes 
invahd. Thereby we enter the general regime discussed in [1], where the break time tbr 
scales with N as 

AtbrxClogAr, (47) 

where constant C depends on the perturbation type. Therefore everything is still 
consistent with the QCC principle. 




5. Conclusions 

In this paper we investigate the correspondence between the classical and quantum 
dynamics of the perturbed cat map on the torus in the hmit of semiclassical small 
perturbations. The correspondence is measured by the overlap between the classical 
density and the Wigner function called quantum-classical fidelity (QCF) and denoted 
by F{t). We study the time evolution of QCF, which stays for a long time at the 
initial value F{t) ^ 1 and than decays towards the ergodic value F{t) ^ faster 
than generally expected. The length of the initial plateau tbr scales with perturbation 
e and Hilbert space dimension as Atbr ~ — log(A^2e), where A is the maximal Lyapunov 
exponent. At the first moment the scaling with seem to be in contradiction with the 
correspondence principle, but this is not the case because the result is only meaningful 
for eN <^ 1. In this particular perturbation regime, the observed behaviour is clearly 
far from general and hence the results presented here for the important and historical 
model of (perturbed) cat map supplement the general knowledge of QCC in evolving 
chaotic systems discussed in [T]. It is important to notice that the presented results can 
be applied to arbitrary chaotic systems which are almost Egorov exact or such that the 
Egorov exactness can be broken by a weak perturbation. 
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